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Abstract 

In this work, a one-dimensional model of crystalline solids based on the Dirac comb 
limit of the Kronig-Penney model is considered. From the wave functions of the 
valence electrons, we calculate a statistical measure of complexity and the Fisher- 
Shannon information for the lower energy electronic bands appearing in the system. 
All these magnitudes present an extremal value for the case of solids having half- 
filled bands, a configuration where in general a high conductivity is attained in real 
solids, such as it happens with the monovalent metals. 
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The application of information theory measures to quantum systems is a sub- 
ject of interest in the last years p!|2|3] . Some relevant properties of the hi- 
erarchical organization of atoms [11[5] and nuclei [6] are revealed when these 
indicators are calculated on these many-body systems. On the one hand, they 
display an increasing trend with the number of particles, electrons or nucle- 
ons. On the other hand, they take extremal values on the closure of shells. 
Moreover, in the case of nuclei, the trace of magic numbers is displayed by 
these entropic products such as the Fisher-Shannon information [7] and a sta- 
tistical complexity measure |8j . Also these statistical quantifiers have revealed 
a connection with physical measures, such as the ionization potential and the 
static dipole polarizability [3] in atomic physics. All of them, theoretical and 
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physical magnitudes, are capable of unveiling the shell structure of atoms, 
specifically the closure of shells in the noble gases. 



A strategy to calculate these indicators is to quantify the discrete hierarchical 
organization of these multiparticle systems through the fractional occupation 
probabilities. These probabilities capture the filling of the shell structure of 
these systems. From them, the different statistical magnitudes are derived. The 
metallic clusters is another system that has also been studied with this method 
jlOj . As in the case of atoms and nuclei, the shell structure of the valence 
electrons is well displayed by the spiky behavior of the statistical complexity 
and the magic numbers are unveiled by relevant peaks of the Fisher-Shannon 
information. 

A different strategy to compute these entropic magnitudes is to use the proba- 
bility density of the quantum system as the basic ingredient. This can be ana- 
lytically obtained in some cases such as the H-atom [3] or numerically derived 
in other cases from a Hartree-Fock scheme [TT][T2] or a density functional- 
theory for atoms and molecules [T3] . 

In this work, we address the problem to calculate these statistical indicators 
in a solid by this last strategy. For this purpose, the band structure of the 
solid has to be determined. The Kronig-Penney (KP) model |T3] is a one- 
dimensional model of crystalline solids that presents a band structure sharing 
many properties with band structures of more sophisticated models. Moreover, 
it also has the advantage that allows to analytically find such electronic band 
structure. 

The KP model considers that electrons move in an infinite one-dimensional 
crystal where the positive ions are located at positions x = na/2 with n = 
±1, ±2, . . ., generating a periodic potential of period a. A simplified version of 
the KP model is obtained when this periodic potential is taken with the form 
of the Dirac comb |T5] : 



where H is the Planck constant, m is the electronic mass and is the intensity 
of the potential. In this case, the spatial part \1'(2;) of the electronic wave 
function is determined from the time independent Schrodinger equation: 



where the eigenvalue E is the energy of the eigenstate "^{x). 

For a periodic potential, the Bloch's theorem [TB] establishes the form of the 
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general solution of the Eq. ([2]). This is a plane wave, with wave number K, 
modulated by a periodic function, uk{x): 

^(x) = e'^^'uKix), (3) 

where uk{x) has the periodicity of the crystal lattice: uk{x) = uk{x + a). It 
implies the following translation property, 

^^{x + a) = e'^^^^x). (4) 

Let us observe that we have the case of free electrons when V{x) = and the 
solutions of type ([3]) recover the plane wave form with a total wave number 
k = K, with uk{x) = constant. It suggests that the solution of Eq. ([2]) in the 
region < a; < a, where V{x) = 0, can be associated in some way with a wave 
number k and then written in the general form: 

"^Ix) = Asinkx + B coskx, (5) 
and, by the translation property (jl]), this solution in the region a < x < 2a is 
^(x) = e'^^lAsinkix - a) + B cosk{x - a)], (6) 
with A, B complex constants and k = ^2mE/h?. 

Two boundary conditions must be fulfilled by ^ at the point a; = a: on the one 
hand, the continuity of the wave function and, on the other hand, the jump 
in the derivative provoked by the Delta function ([1]). This gives the relations: 



^(a + 0) = ^(a-0), (7) 
^'(a + 0) = ^'(a-0) + 2fi^(a). (8) 

From these boundary conditions applied to the wave functions ([MS]), the fol- 



lowing homogeneous linear system is obtained for the unknonws A and B: 



( 



sin ka 



cos ka — e*^"" 



[ke — k cos ka — 20 sin ka) [k sin ka — 2Vl cos ka) 




■ (9) 



To have a non-trivial solution, the determinant of this 2x2 matrix has to be 
zero. Then, the following quantization relation for k is obtained jl5j : 



cos Ka = cos ka -\ sin ka . 

k 



(10) 



The electronic band structure of the one-dimensional crystal is contained in 
this equation. When K varies its value in the different Brillouin zones, given 
by (m — l)7r < \Ka\ < rnvr, with m = 1 for the first Brillouin zone, m = 2 for 
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the second Brillouin zone, etc., only certain intervals of k are allowed. These 
intervals for k are the energy bands of the electronic system. The positive and 
negative parts of these intervals correspond with the positive and negative 
parts of the Brillouin zones, respectively. In the limit Qa = 0, the free electron 
problem is recovered, then the solutions are the plane waves with k = K. In 
the limit Qa = oo, we have the square well problem, then the wave number of 
the eigenstates verify sin ka = 0. For an intermediate case, < Qa < oo, Eq. 
f lTOj) has to be solved. Concretely, for the particular value fla = 4, that has 
also been used in Ref. [TS], the lower energy bands obtained in this system for 
k > are: 



The bands are symmetrically found for k < 0. Observe that, to finally get the 
wave function of the electronic states, we additionally need the normalization 
condition to completely determine except a global phase factor. For our 

calculations, by taking Qa = 4, we will perform this normalization in the unit 
cell [0, a]. 

The basic ingredient to calculate the statistical entropic measures in which we 
are interested is the probability density of the electronic states. This is given 
by p{x) = From this density, we proceed to compute the statistical 

complexity and the Fisher-Shannon information. Notice that the wave function 
for a given k is transformed in — for —k, therefore all the magnitudes 
depending on the density are the same in both cases, and then we reduce our 
study to the positive part {k > 0) of the electronic bands. 

Let us recall the definition of the statistical complexity C piTT] , the so-called 
LMC complexity, that is defined as 



where H represents the information content of the system and D gives an idea 
of how much sharp is its spatial distribution. As a quantifier of if, we take a 
version used in Ref. ^7j. This is the simple exponential Shannon entropy [18] , 
that has the form. 



2.154 <ka< TT 
4.578 <ka<27T 
7.287 <ka<37T 
10.174 < ka<ATT 



(1st band), 
(2nd band), 
(3rd band), 
(4th band). 
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C = H -D 



(12) 
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Fig. 1. Shannon entropy, S, vs. the adimensional wave number, ka, for k > 0. Only 
the four lower electronic bands given in expression (jlip are shown. 

For the disequilibrium we take the form originally introduced in Refs. ISfTTj . 

that is. 



The entropy, S, and the statistical complexity, C, for the lower electronic 
bands of the present one-dimensional crystalline solid are given in Fig. [1] and 
Fig. [2], respectively. When this hypothetical solid is in a situation of high 
conductivity, i.e. when it contains a band that is partially filled and partially 
empty, it can be observed in the figures that the more energetic electrons 
attain the highest entropy and the lowest complexity in the vicinity of the 
half-filled band. This is the point where in general the highest conductivity is 
also attained. Take, for instance, the real case of the monovalent metals, that 
include the alkali metals (Li, Na, K, Rb, Cs) and the noble metals (Cu, Ag, 
Au). These metals present all the bands completely filled or empty, except 
an only half-filled conduction band [20]. Compared with other solids, these 
metals display a very high conductivity, that in the cases of Ag and Cu it 
is the highest in nature. Then, it is remarkable this coincidence at the point 
of half-filled band where, on the one hand, the entropy and the statistical 
complexity are extrema for this model of solids and, on the other hand, the 
conductivity reaches its upper values for the real cases of monovalent metals. 

Now, we check that other statistical entropic measures also display this be- 
havior when the solid has half-filled bands. Let us take, for instance, the 
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Fig. 2. Statistical complexity, C, vs. the adimensional wave number, ka, for A; > 0. 
Only the four lower electronic bands given in expression (jlip are shown. 

Fisher- Shannon information, P, that has been apphed in different contexts 
[rpT] for atomic systems. This quantity is given by 



J-/, 



where the first factor 



J 
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(17) 



is a version of the exponential Shannon entropy [18], and the second factor 

f« [dp{x)/dx]'^ 



dx 



is the so-called Fisher information measure |22], that quantifies the stiffness 
of the probability density. Observe in Fig. [3] the confirmation of the previous 
results obtained in Figs. [1] and |2] for S and C, in the sense that the extremal 
values of P for this model of solids are also reached at the half-filling band 
points. 

The former calculations are done orbital by orbital, i.e. thinking that the solid 
is a set of individual and independent orbitals, each one identified by its own 
wave number k. We can change the point of view of the problem and to think 
that the solid stands in some kind of collective state whose probability density 
Pt{x) is the normalized sum of all the allowed electronic densities obtained from 
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Fig. 3. Fisher-Shannon entropy, P, vs. the adimensional wave number, A;a, for /c > 0. 
Only the four lower electronic bands given in expression (llip are shown. 

the orbitals with wave numbers in the interval [kmin,kmax]] ^mm will be the 
minimal electronic wave number of the solid, i.e. the lowest k obtained in the 
first band, and kmax will be the upper k corresponding to the most energetic 
electron of the solid. The expression for pt{x) is 



Observe that pt{x) is normalized in the interval [0, a], Jq pt{x)dx = 1, and that 
in the present model of solid kmin = 2.154 as given in formulas (fTTi) . 

The calculation of the statistical complexity Ct for this pt{x) is presented in 
Fig. m In this case, the minimal values of Ct are also located in the vicinity 
of the half-filled electronic bands such as the behavior of C for the individual 
orbitals shown in Fig.O In the hypothetical limit case of a solid where kmax^ ^ 
1, let us remark that the density pt{x) will tend to the uniform density and 
then the lowest value of complexity, Ct = 1, can be reached, as it can be seen 
in Fig. m 

In summary, this work puts in evidence that certain conformational prop- 
erties of many-body systems are refiected by the behavior of the statistical 
complexity C and the Fisher-Shannon information P. In the present study, 
the electronic band structure of a model of solids has been unfolded and the 
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Fig. 4. Accumulated statistical complexity, Cj, vs. the adimensional wave number, 
ka, for A; > 0. Only the four lower electronic bands given in expression are 
shown. 



measurement of these magnitudes for such a model has been achieved. It is 
remarkable the fact that the extremal values of C and P are attained on the 
configurations with half-filled bands, which is also the electronic band config- 
uration displayed by the solids with the highest conductivity, let us say the 
monovalent metals. Therefore, the calculation of these statistical indicators 
for a model of solids has unveiled certain physical properties of these systems, 
in the same way that these entropic measures also reveal some conformational 
aspects of other quantum many-body systems. 
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